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ABSTRACT 


The Airborne Digital Data Acquisition System (ADDAS) was developed by the 
joint Naval Postgraduate School/Harvard Unversity research project, called the Ocean 
Prediction Through Observation, Modeling, and Analysis (OPTOMA) Program, for the 
real-time synoptic mapping of oceanic features. 

Given the capabilities of such a svstem. an aerial sampling strategy was 
developed based on the spatial scales of oceanic variables and the aircraft and ADDAS 
limitations. The strategy was applied to “simulated” AXNBT survey flights on 
“synthetic” fields. The simulated AXBT observations were cbjectively analyzed and 
compared co the synthetic fieid. By varying the simulated AXBT spacing and the 
objective anaivsis parameters, an “optimum” AXBY spacing to adequately representma. 
synthetic structure of the domain was deduced. 

To further develop such a svstem ior both airborne and shipboard use, a 
proposed Shiv/Aircratt Data Acquisition. Dispiay, and Analvsis System (SADADAS) 
was oresented. ‘the tollowing synoptic mappings produced could provide a possibie 
tactical decision aid for antisubmarine warfare (ASW) operations: 1) Sea surface 
temperature; 2) Temperature/Sound speed at selected depths; 3) Maps of mixed layer 
depth; 4) Subsurface ducts/sound channels; and 5) Vertical sections of 


temperature/sound speed. 
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I. INTRODUCTION 


The development of an Ocean Descriptive/Predictive System (ODPS) for the 
study and forecasting of oceanic synoptic mesoscale features and the study of the 
California Current System are the goals of the Ocean Prediction Through Observation, 
Modeling, and Analysis (OPTOMA) Program, a cooperative effort between the Naval 
Postgraduate Schooi (NPS) and Harvard Universitv. The Airborne Digital Data 
Acquisition System (ADDAS) was deveioped at NPS in response to the need for a 
real-time, large-scale synoptic mapping of the mesoscale features for both the 
initialization of dynamic model forecasts and the verification of the dynamic models 


used. 


pee OBJECTIVES 

To deveiop a system sucn as the ADDAS for use in operational antisubmarine 
warfare (ASW) is the primary objective of this thesis. With the capabilities of such a 
system: 1) What optimal aeriai sampiing strategy should be utuized to acquire oceanic 
data for naval planning and operations? 2) What further enhancements could be 
developed for a system like the ADDAS, and how could such a system be applied to 
shipboard operations as well? 3) What operationally applicable synoptic fields can be 
developed from the data, and how could such fields be applied in planning and 


execution of Naval operations, specifically in ASW? 


B. BACKGROUND 

With the availability of expendable aircraft deployable instruments such as 
airborne expendable bathythermographs (AXBTs), sound speed profilers (AXSVPs), 
current profilers (AXCPs), and soon-to-be-available air-launched conductivity, 
temperature, depth profilers (AXCTDs), the aircraft is the platform of choice for near 
real-time data acquisition. The P-3 Orion ASW patrol aircraft supports a system such 
as the ADDAS exceptionally well. It has the payload capacity and avionics for air 
launched sensors as well as substantial on-station endurance and all-weather capability. 
It can cover a sizable oceanic region in a matter of hours as compared to days for a 
research vessel. The ADDAS has been deployed on six OPTOMA P-3 missions to date 


and is a proven (research) operational system. It has been invaluable as a primary 
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means of data acquisition for the OPTOMA Program. Additionally, carrier-based 
ASW aircraft such as the S-3 Viking, and rotary-wing ASW aircraft including the SH-2 
Sea Sprite, SH-3 Sea King, and the SH-60 Seahawk, can provide additional airborne 


sensor deployment platforms. 


C. MOTIVATION 

The capability for near real-time mapping of oceanic variables has a direct 
application to ASW operations. As submarines become quieter, the warfare 
commander will have an even more difficult task of localizing and tracking them. An 
increased knowledge of the submariners’ environment ‘vill be required. An AXNBT 
survey flight can provide a “nowcast” of the region for use in tactical briefings, and 
such planning should be a key factor in ASW operations. With maps from such 
AABT surveys, the etfects of the ocean structure On standard tactics can Gemiddne sam 
and appropriate modifications ‘o tactical doctrine can be suggested. Additionally, 
subsequent survey tlights could be compieted to identify the temporal and spatial 
changes over the course of the exercise or operation and to retresh the original Sup vere 
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‘Traditionally, ASW tactics and training are based on a “generic” ocean structure. 
and usually only one or two AXBT’s are deploved in an exercise which tendsmee 
reenforce this concept. One may argue that a single AXBT could be representative of 
the ocean thermal structure in some region. However, if an AXBT was deployed in the 
middle of an eddy or in an area of strong temperature gradients, as are observed in the 
Gulf Stream, it may not be representative of that region. In the research field, the 
ADDAS has proved itself as a powerful tool. It is time for this capability to be 
developed for operational use in the fleet. 

The current version of the ADDAS and similar systems in use today are 
presented, followed by the oceanographic application of the ADDAS in the OPTOMA 
Program. Next, an optimal aerial sampling strategy is developed and applied to 
“synthetic” oceanic domains and compared through use of the objective analysis 
technique. A proposed Ship/Aircraft Data Acquisition, Display, and Analysis System 
(SADADAS) concludes, with application to both land-based and shipboard aircraft 


operations. 
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Il. THE AIRBORNE DIGITAL DATA ACQUISITION SYSTEM (ADDAS) 


A. ADDAS DEVELOPMENT 

The OPTOMA Program’s Airborne Thermal Structure Mapping Project at NPS 
culminated with the development of the ADDAS and its successfui research use. [In 
the early survev. flights of this project, AXBT analog traces were recorded onooara the 
P-3 aircraft by the AQA-7 iofargram svstem recorder. The iraces were nand digitized 
using gridded templates based on the standard Navy frequency-to-iemperature and 
fail-rate-to-depth equations. This produced traces of coarse resolution and resulted in 
a substantial time-lag for the generation of “real-time” synoptic maps. 

eer ew pemoes  .nstilute Oi Occanograpny, Vir. Meredith H. Sessions iad 
Meweloped a@M@sVvStem .:or <he recording ana real-time digitai conversion of an AAST 
analog signal. which nad oeen used successtuily in the North Pacific Experiment 
(NORPAA) (Sessions, 1980). The Scripps system was made availabie to the OPTOMA 
feoerdm) and was used :or the cost-flight digitization o! unalog tape reccrdings of the 
ANDI signals acquired by the OPTOMALIP ‘hght on i8 Julv (984. inspired oy the 
relative ease Of data acquisition and the improved resolution of AXBT profiles 
provided by the Scripps system, as well as the recent development of an AXBT 
digitizing circuit for the Sippican MK9 digitizing unit, a prototype ADDAS became 
feasible and was developed at NPS (Colton and Mooers, 1985). First deployed on the 
OPTOMAI3P flight on 27 October 1984, the ADDAS has proven invaluable in 


acquiring real-time data for synoptic mappings and model initialization. 


B. THE PRESENT ADDAS 

The ADDAS is a self-contained, semi-portable data acquisition, processing, and 
analysis system. It utilizes off-the-shelf hardware and modified shipboard XBT and 
CTD data acquisition and analysis software developed by the OPTOMA Program at 
NPS. The heart of the system is a HP9816 MC68000 16-bit microprocessor with a 
HP9121 dual 3.5 inch microfloppy disk drive. Aircraft 28 VDC power is inverted to 
115 VAC, 60 Hz power by an internal static inverter. AXBT drop positions can be 
entered either manually by the operator through the computer keyboard (latitude and 
longitude, as determined by the Tactical Coordinator/Navigator from aircraft inertial 
navigation) or automatically by the ADDAS’s optional LORAN C unit. The AXBT 
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audio signal is received through the ARR-72 aircraft sonobuoy receiver, digitized in 
real-time by a Sippican MK9 front-end processor, and the trace is displayed on the 
computer monitor. Once the AXBT has completed its fall, the digitized profile is 
stored on a data diskette. Currently, the ADDAS is capable of digitizing only one 
AXBT at a tume. To simulate a multi-channel capability, one channel of a stereo hi-fi 
video cassette recorder (VCR) is used to record the simultaneous audio signal from a 
second AXBT (of different VHF carrier channel) for post-flight digitization and 
analysis. The audio signal of the digitized AXBT is also recorded on the other VCR 
channel as a backup. An oscilloscope provides a visual display of the AXBT audio 
signal to nelp detect faliures both in the aircrait sonobuov receiving svstem and in the 
AXBTs. Once the survey is complete, maps of objectively anaivzed variabieS are 
oroduced on the monitor and hard-copy output is generated ov a HP7470A two-pen 
plotter, either on the return flight or on the ground. A list of ADDAS aardware 


somponentsand zheir @urrent retail COSt 1s Includea@ mie pend 


C. ADDAS LIMYTATIONS 


Limutations to the svstem can be classed into two Categories, iniiisicueaes 


sxtrinsic. These limitations affect ADDAS operation and ultimately the Gensityaiem 
ACN Bd SOVERICe Aa cuey, 


|. Intrinsic Limitations 
These limitations are specifically inherent in the system design and operation. 
They are: 

e Not multi-channel capable (presently). Two AXBT’s of different channels cannot 
be digitized simultaneously. An additional MK9 or the development of a 
multi-channel MK9 would be required. In addition, an advanced microprocessor 
wae eS Ue capability is required. These are very feasible with todays 
technology. 


e No hardware redundancy. Although the system has. proven exceptionally 
reliable, if one component should fail the real-time capability will be lost. 


e Inverter operation. The inverter is rated at 5.6 amps. With only the 
microprocessor, disk drives, MK9, and VCR in the circuit, the inverter operates 
at 6.3 amps continuous. This has not, been,a problem to date since the VCR 1s 
only used during data acquisition and it 1s disconnected when the plotter is being 
used. An inverter with a higher current rating may be required. 

2. Extrinsic Limitations 
These limitations affect the overall capability of the system to perform its 
given mission and are related to the sampling platform and the sensors employed. 
Tineveacc: 


e Aircraft weight and balance. Aircraft are constrained By a maximum gross 
weight for take off. It is a function of the environmental conditions at the 
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departure point. For navy aircraft, these limitations are specified in the 
respective NATOPS oyava) Air Training and Operations Procedures 
Standardization) Manual for each aircraft type and govern safety of flight. In 
order to carry enough fuel to complete the mission, and have the required reserve 
fuel at ee there may need to be a reduction in the number of AXBT’s or 
other sensors fo be deployed in order to remain within NATOPS limitations. 


e Aircraft airspeed and altitude. These factors affect aircraft fuel consumption. 
Aircraft engines operate_more efficiently at high altitude yet the ADDAS requires 
a lower altitude for effective operation, A compromise is made. | Nr. 
transit to and return from the survey region will 1s accomplished at high altitude 

> 10,000 ft MSL) and at an optimum speed (~310 Kts indicated airspeed). 

nce on-station, the aircraft will be operated at loiter airspeed (180 to 210 kts) 
and altitude will normallv be between 500 and 1000 ft AGL (above ground leve!). 
This will cover as large an area as possible, and provide for optimum ADDAS 
performance. Ajircratt aititude affects the time tt. takes tor the AXBT. to tmpact 
the water after deployment from the aircraft and the line-or-sight radio 
transmission range. At nigher altitudes. the transmission range is_incréased and 
interference from other sensors on the same RF channel may aifect the AXBT 
transmission. Airspeed and altitude must be considered in fuel planning (weight 
and balance) for the mission. Additionaily, safety of flight requirements as 
specified by NATOPS and squadron SOP (standard operating procedures) may 
require higher altitudes for certain evolutions such as pilot changes. 


® AXBT fall duration through che water; for shallow (300 m) and deep (760 m) 
Post? Saee0O and S00 seconds. respectively (Sessions, ef ai.. i976, Bane and 
Sessions, 1984). | 

¢ Number of AX3T channeis. Oniy three channels, {2. !4, and 16 (corresponding 
foe 2 Oana.) Vilnz, respectively.) are currently manufactured «Gent, 
1982}. Two or more AXBT’s of the same channel_musi not oe deploved ‘vithin 
the line-of-sight Wr transmission range. [he contlicting signals can trigger tne 
Vin? and rcsuit m the 10ss of a prolile. it is basicaily a tunction of aircralt 
milituderand speed. 


Mexccinal radio imterrerence on AXBIT carrier channels. Once, the M9 is 
initialized and awaiting the AXBT signal, spurious radio transmissions on the 
receiving channel can pre-trigger the MK9 and the digitized profile may be lost. 


D. SIMILAR SYSTEMS IN USE TODAY 

In addition to the Scripps system, the Naval Ocean Research and Development 
Activity (NORDA) also has a similar system, the Airborne Data Acquisition and 
Processing System (ADAPS) (Miles, 1984). ADAPS is currently used by 
NAVOCEANO for work in acoustics. These systems were designed for specific 
missions and are multi-channel capable. The organizations have the advantage of a 
technical staff and maintenance support, and as a result these systems are 
one-of-a-kind in design and operation. The NPS ADDAS, in comparison, utilizes 
totally off-the-shelf hardware. It also has the added capability of being employed in 
shipboard XBT surveys as well with a modified software package and different 
digitizing circuit board for the MK9. However, the ADDAS is limited in that it is only 
single-channel capable at the present. With further advances in expendable sensors, 


off-the-shelf digitizing units, and computer systems capable of multi-tasking, systems 
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such as the ADDAS, ADAPS, and the Scripps system can become available to the 
Naval Oceanographer and research community to more fully describe the oceanic 


environment. 
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Hil. OCEANOGRAPHIC APPLICATION OF THE ADDAS IN THE 
OPTOMA PROGRAM 


Near-real-time synoptic mappings of the mesoscale eddy field are a Key facet in 
the observational network of an ODPS. These mappings, using data acquired through 
the OPTOMA Program’s AXBT survey flights, provide nowcasts of the mesoscale 
variablity in the OPTOMA domains which are used for shipboard survey dlanning, 
dynamic model initialization. modei forecast verification. and assessment of the 
interannual variability in the domains. 

The NOCAL (Northern California) and CENCAL (Central California) 


— 


Bao@omains are the major survey regions for the OPTOMA program (figure 3.1). 
eseuesUodomain 1s an open ocean region in the Caufornia Current Svstem (CCS) 
meer 'l0-km Square ‘with water depths .f the order of 4+ km. ihe NOCAL 
subdomain, iS0-km downstream from the Mendocino Escarpment and about iS0-km 
Offshore, is the primary survey subdomain. [ts bottom topography is characterized ‘ov 


gently sloping seatloor and no seamounts, iithougn seamounts, the Mendccino 
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Eeeanpment, and the continental margin are nearbv. NOCAL is the simpiest domain 
that can be found off California for a region this size (Mooers, 1986). AXBT survey 
flights generally expand the subdomains to about 260-km square in order to acquire 


sufficient data to initialize the boundary conditions for model runs. 


A. AXBT SURVEY FLIGHT PLANNING 

As with any oceanographic survey, planning is essential for the successful 
completion of the mission. Formal tasking for research flights is requested through the 
chain of command from the research sponsor. AXBT survey flights are then scheduled 
on a not-to-interfere basis and are contingent upon the operational tasking of the 
patrol wing. Once flight time has been authorized, coordination with the Patrol Wing 
Schedules Officer, preferably two to three months in advance of the date of the 
proposed survey flight, should begin. Flexibility is the key, and periodic contact should 
be maintained at least bi-monthly as changes in operational commitments (within both 
parties) may require the rescheduling of the flight. One month prior to the flight, a 
draft flight plan is prepared, and required mission personnel, supplies, equipment, and 
transportation are identified. Within two weeks of the flight, the Schedules Officer 
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subdomains. 
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should know which squadron will be tasked to provide the services. Coordination is 
now primarily with the squadron. 

The proposed flight plan should be discussed with the squadron Operations 
Officer and Schedules Officer. Changes for safety-of-flight, number of AXBT’s to be 
deployed, and the planned flight itinerary should be included in the revised flight plan. 
Additionally, squadron personnel should be requested for unloading and preparing the 
AXBT’s for the flight and loading them onto the aircraft. 


B. AXBT SURVEY FLIGHT EXECUTION 

OPTOMA survey flights nave generally deploved between 70 and 96 AXB'T’s per 
mission. Once on station in the survey subdomain, AABT’s are normally launched at 
15 nautical mile (NM) intervais, alternating shallow and deep AXBT’s. The SSQ-36A 
AXBT is manufactured in two factory-set depths. 505 m (currently the Navy standard 
AXBT) and 760 m (pnmarily used for researcn). The shallow AXBT’s are digitized in 
real-time and stored by the ADDAS, while the deep AXBT’s are recorded on tape ior 
post-flight digitizing. A flight log is maintained with specific information about each 
Grop: <ime, latitude, longitude.remarks. All shallow AXBT’s are of one channet, while 
mime cesp ANB s are of a different channel. This limits the radio interference 
memvesen AAG! drops since only two channeis. i4 and |6, can woe used in the 
OPTOMA survey regions. (Channel 12 is not used as it 1s the California Forest 
Service emergency broadcast channel.) 

With the survey completed, the stored AXBT profiles are read and the desired 
output variables; e.g., sea surface temperature (SST), mixed layer depth (MLD), depth 
of the 8°C isotherm (D8C), sound speed at a chosen depth (CSZ), or dynamic height 
relative to a specific depth (DHZ), are selected. These variables are then entered into a 
statistical objective analysis routine based on an application by Carter and Robinson 
(1981). A discussion of the objective analysis technique as it is applied in the ADDAS 
software is presented in Chapter IV. This technique uses a two-dimensional correlation 
function based on the statistics of the data set. The five most correlated values within 
a 55 km radius of an interpolation point are then used to compute a value at the 
interpolation point. Contoured maps of the selected fields are then generated and 
printed on the plotter. These plots are available prior to landing so that an immediate 
debrief can be accomplished. Upon return to NPS, the recorded deep AXBT’s are 
digitized and added to the data set. The objectively analyzed fields are then generated 
on the IBM 3033 mainframe using both the shallow and deep AXBT’s. 
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The most recent AXBT survey flights were OPTOMAI8P consisting of two 
consecutive flights on 31 October and 2 November 1985, and OPTOMA20P on 16 
March 1986. Summaries of each mission and the data acquired are presented below. 

1, OPTOMAI8P 

Aircraft and flight crews were provided by Patrol Squadron Forty-Six, 
assigned to Commander, Patrol Wing Ten, for the 31 October flight and Patrol 
Squadron Ninety-One, assigned to Commander, Reserve Patrol Wing Pacific Fleet, for 
the 2 November flight. Both are based at NAS Moffett Field. The Sight tracks were 
executed as planned (Figure 3.2). A total of 180 AXBT’s were depioved on the two 
flights. Spacing was at 15 NM intervals, alternating deep (Ch. 16) and shallow (Ch. 
14) AXBT’s. On-station ground speed averaged 200 kts for each flignt and the altitude 
varied between 500 and 800 ft AGL, depending on the low-level visibility. The logged 
positions were obtained from the aircraft inertial navigation system with an accuracy of 
PVG Kaa). 

d. iMission Critique 

Two different models of the P-3 aircraft were utiuzed on the flights: “aaa 
on 31 October. and a P-3B on 2 November. The following differences are igiam 
oetween the two models tor \DDAS installation: 

¢ The 28VDC power connections in the P-3C were made to terminal board 426 
ae eee aft of sensor station 2 (SS2) bulkhead. The terminals used 
were CIO (+28VDC_ power_SS2) and Cll (ground). In the P-3B, power 
connections were made to TB431 located above SS2.° Terminal, A24. provided 
+28VDC power while the ground was a stud located on the far right side of the 
terminal board. The sonobuoy pin connections were the same for both models of 
aircraft. (Please see Colton and Mooers, 1985, for a complete description of 
ADDAS installation.) 

¢ The P-3C model does not have a LORAWN system incorporated into its avionics 
package. Thus, no LORAN antenna was available to be connected to the 
ADDAS LORAN unit, so all AABT Bae were obtained from the inertial 
navigation system and checked against the OMEGA radio navigation unit. 

e The physical internal aon of the P-3C is different from the P-3B. 9 diierems 
much less room for the ADDAS and operator in the P-3C at sensor station 2 due 
to the bulkhead separating stations. e ADDAS operator had to work literally 
“shoulder to shoulder” with the sensor station | operator. 

The 31 October flight with Patrol Squadron Forty-Six went very well. 
Since the ADDAS LORAN was not used, the positions of each AXBT drop were 
taken from the inertial navigation and entered manually into the AXBT program. 

The 2 November flight proceeded very smoothly due to the experience 
gained by Patrol Squadron Ninety-One from previous OPTOMA flights. ADDAS 


installation was completed in a record time of forty minutes. AXBT positions were 
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Figure 3.2. Flight Plan, OPTOMAISP. 
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obtained as on 31 October due to problems encountered on previous flights with RF 
interference from aircraft radios affecting the LORAWN signal. The only notable 
difficulty on this flight was that a low-level (500 ft) stratus deck required an on station 
altitude of 800 ft AGL for flight safety. This increased the transmission range of the 
AXBT’s. On the first southbound leg, a tailwind increased aircraft groundspeed to 235 
kts and resulted in pre-triggering of the MK9 due to overlapping AXBT signals on the 
same channel. On subsequent downwind legs, airspeed was adjusted to maintain 200 


kts groundspeed. Table | lists a summary of AXBT expenditures for the two flights. 


TABLE 1 
AXBT EXPENDITURE SUOMI Via i te 





rhight Area Cuanne! Denioyed Caiied | Recovered | 
SL OCe “SEC 6 = ) a . 
.~ Nov NOC Is 6 46 6 <0 | 

5 =5 - 4} ! 
Two-day Totai: 180) a 160 


Notes: 


' Due to no radio signal transmission and iate iight-olf of the probe. 





b. AXBT Survey Results 

The uncorrected positions of the shallow AXBT’s, from the 31 October and 
the 2 November flights, and the analyses of D8C, SST, and CS150 that were generated 
onboard the aircraft on the return leg from the 2 November flight were examined 
(Figures 3.3a, 3.3b, 3.3c, 3.3d, respectively). A compiled version of the objective 
analysis program, using a floating-point accelerator was applied. Only the shallow 
AXBT data (30 NM spacing) were combined to produce the analyses. The data were 
used as received and digitized by the ADDAS without editing or filtering performed on 
profiles in order to provide a real-time analysis on the return flight to NAS Moffett 
Field. Some quality checks are made by the ADDAS software. A complete summary 
of the acquired AXBT data is presented by Wittmann et al. (1986). 
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Figure 3.3d Sound Speed at 150 m (m/s), as in igetnesea. 
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The inflight digitized channel 14 AXBT data were later combined with 
post-flight digitized channel 16 AXBT at NPS. These data were stored and processed 
on the IBM 3033 mainframe computer. Profiles were edited to remove spikes, and bad 
profiles were eliminated. Post-flight, the recovered AXBT locations were examined 
(Figure 3.4a). 

¢ D8C, meters (Figure 3.4b): This analysis presents a “topographic” map of the 
depth of the 8°C surface, which is genérally accepted as a good indication of tne 
location of the permanent thermocline and pvcnocline for this region. From the 
D&C analyses, the streamfunction for the surface geostrophic flow can be inferred 
with good accuracy. The major features were two cyciones and two anticyclones 
(indicated by “Cs and “A’s’, respectively), alternating along the coast. Ine 
cycione at the far north and the anticyclone at_the far south were fairly weak 
compared to the cyclone ciose to the coast near Point Reyes (diameter about 1U00 
km, depth variation of over 60 m) and the stronger anticyclone off Point Arena. 
(diameter, about [50 km, depth variation of over 60. m). The juxtaposition ot 
these eddies contributed to alternating onshore and offshore flow. 


@ SST, °C (Figure 3.4c): The main features are two cold core centers findicated ov 


“s S ) near the coast and a generailv weak positive offshore temperature gradient. 
The cold core center in the south or the sampled region was elongated, parailei to 
meee COdsc Wilercas “ilat im) ihe morth of the sampled region ‘near Cape 
Mendocino) was more circular. stightlv elongated offsnore. 

* CS150, m/s (Figure 3.4d): Sound speed was computed using a constant salinity 
Of 53.9 Dp. The sound speed '!s fairly constant with oniy a 6 ms varianiuity, 
aithough there is an obvious_correlation oetween sound speed patterns and the 
pocn oatterns 0! the Gepth of the $ degree isotherm. the strong feature evident 
ier tsure 5.50 JUSt tO the south of the northernmost cold center was a orobiem 


a 


with the data which was eliminated in the post-flight data processing. 
2, OPTOMA20P 

A combined flight crew from the Naval Air Reserve Center and Patrol 
Squadron Ninety-One at NAS Moffett Field performed the mission. The flight plan 
was executed as shown in Figure 3.5. The ten flight tracks were spaced 15 NM apart. 
Ten AXBT’s were deployed along each track at a 15 NM spacing, with the exception 
of only eight AXBT’s deployed on tracks 1 and 10 (no AXBT’s were deployed at the 
corners of the survey region), for a total of 96 AXBT’s. Groundspeed averaged 220 kts 
and on-station altitude varied between 500 and 2000 ft AGL for safety of flight due to 
turbulence and scattered convective activity and thunderstorms in the area. 

AXBT’s were to be deployed alternating channel 16 (deep, 760 m probe) and 
channel 14 (shallow, 305 m probe). Due to inventory error, all AXBT’s for this flight, 
both channel 14 and 16, were shallow (305 m). This was not realized until the first few 
channel 16 AXBT’s were deployed and their time of fall observed. AXBT’s were then 
launched at approximately four minute intervals, alternating between channel 16 and 
14 AXBYT’s. AXBT positions were obtained from the aircraft inertial navigation system 


(2 NM accuracy) and entered in the flight data log. The positions were not corrected 
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Figure 3.4a Recovered AXBT Drop Locations, OPTOMAIS8P. 
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for navigation error. The channel 14 positions were manually entered in the AXBT 
acquisition program and the signal was digitized in real time by the ADDAS and 
stored on diskette. Both channel 14 and 16 AXBT analog audio signals were recorded 
On magnetic tape using a stereo hi-fi VCR. Channel 16 AXBT’s were digitized on 
return to NPS. 

c. Mission Critique 

The flight proceeded very well even though the weather was less than 
optimum. The survey region ‘vas dominated by iow stratus and cumulus clouds as weil 
as a few towering cumuius. Moderate turbulence was experienced throughout the 
duration of the flight. Additionaliy, scattered rain squalls required deviation from the 
planned flight tracks for safety of fight. Winds at flight level were southeastward at 30 
knots. and the sea state was 5. ‘The navigation team did an excellent job and the 
resulting AXBT drops ended up very nicety spaced and close to planned drop positions. 

A maifunctioning channei 16 sonobuoy receiver on the aircraft distorted 
the AXBT audio carner signal. The ADDAS oscilloscope provided a visual reference 
femetme inilignt Technician (TP 1} to troubleshoot and repair the receiver. fortunatety, 
only 8 channel 16 AXBT audio signals were lost and could not be digitized oy the 
Beepican Mark 9 unit ipon veturn to NPS. Due to a minor problem with ‘he 
lofargram recorder, only 6 of the 8 AXBT traces were obtained and these profiles were 
hand digitized. 

The use of a stereo hi-fi video cassette recorder (VCR) in the audio mode 
proved to be a worthwhile addition to the ADDAS system. With its capability of 
recording up to six hours of audio on a T-120 VHS cassette, the operators only had to 
change the tape once as compared to previous flights where the AXBT audio signal 
was stored on standard audio cassettes with only 45 minutes of recording on each side 
of the tape. No profiles were lost during tape changes or because a tape ran out. The 
VCR was mounted in the ADDAS cabinet on a layer of 3/4-inch foam rubber, which 
helped insulate it from vibration caused by turbulence encountered during the flight. 
No degradation in the recorded signal was noted. 

The internal ADDAS inverter is rated at 5.6 amps. With the inclusion of a 
VCR in the system, the inverter was running at 6.3 amps continuous. This did not 
appear to be a problem and the inverter did not overheat or malfunction. To be on 
the safe side though, once data acquisition was completed, the VCR was disconnected 


from the circuit when the plotter was to be used. 
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The Sippican MK9 digitizes the AXBT audio signal ten times each second. 
MK9 output was sampled for three-tenths of a second, producing three data points. 
The median value of the three digitized points was stored. This step was inciuded in an 
attempt to eliminate non-representative spikes in the temperature trace caused by RF 
or other interference in the signal. This feature appeared to work well, with the 
exception of long duration breaks in the signal, possibly caused by sea wash-over of 
the AXBT transmitter. A minimum of post-flight editing was required on the profiles. 
(A compilation of the AXBT data is presented by Ciandro er al., 1986.) Table 2 lists a 


summary of AXBT expenditures for this mission. 


TABI 
AXBT EXPENDITURE SUMMARY OPTOMA20P 


is 


Channei Denloved fatied* Osis Deteted” Recovered | 
| 
(4 48 t ) 0 AT | 
Lo <8 0 2 | 39 | 
otal: 96 f i 86 
NIQUES: 


| Due to no radio signal transmission and late light-off of the probe. | 


* Profiles lost due to absence of lofargram trace from which to hand digitize. 


3 Bad profiles deleted during post-flight editing. 





b. AXBT Survey Results 

A representative sample of AXBT traces show the variability of the thermal 
structure in the survey region (Figures 3.6a and 3.6b). The AXBT drop positions, and 
the analyzed fields of D8C, SST, and MLD were produced by the ADDAS on the return 
leg of the flight (Figures 3.7a, 3.7b, 3.7c, 3.7d). Only channel 14 AXBT data (30 NM 
spacing), as received and digitized by the ADDAS, were used. The fields were 
generated by a compiled version of the objective analysis program using a 
floating-point accelerator. 

Post-flight digitized channel 16 AXBT data were combined with the inflight 
digitized channel 14 AXBT data (Figure 3.8a). These data were stored and processed 
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Figure 3.6a 
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Figure 3.7a. Shallow AX T Drop Locations, OPTOMA20P. 
Origin 38°N, 125°W. Tick marks = 25 km. 
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Figure 3.8a Recovered AXBT Drop Locations, OPTOMA20P. 
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on the IBM 3033 mainframe computer at NPS. Remaining spikes were removed and 


bad profiles were deleted from the data base. 


e D8C, meters (Figure 3.85): The most notable features are the major ¢Cvclemms 
‘counterclockwise, indicated by a “C” ) center in the northwest with its associated 
ow zone extending to the south-southeast, as indicated bv the shoaler depths, 
and the. very large anticyclonic (clockwise, indicated, by an “A”) center 
characterized by the depression in. the southeast. Additionally. half of two 
smaller cyclonic flows aré_observed in the north as well as a small cyclone near 
the center of the region. These features indicate an onshore jet off Pt. Arena. 


® SST, °C (Figure 3.8c): Major features of interest are the cool zonemmaumes. 
northwest and the warm center (indicated 5v a “W”), corresponding to the large 
evcione inthe northwest and the iarge anticvclone to the southeast, respectively. 
ihere is also evidence of a cool center ("XK ) 1n the north corresponding (o/Omemes 
the two smaller cyclonic tlows. 


* MLD, meters (Figure 3.sd): Of interest in this analysis is the wide vangemem 
values for MLD, irom less than 16 mto more_than 68 m. The mean ViLD 4s 30 
m. This variapility could be, the resuit of the large atmospheric occluded 
loW-pressure system that dominated the region on the previous day. As the 
system passed, the winds shifted trom generaily northeastward to southeastward 
at 20 to 30 knots and the surtace laver was apparently im a State Of rapid eam 
during the survey period. All AXBT profiles showed a very isothermal surface 
layer. There was some association of the patterns in MLD with those of DaG@ 
for exampie. zne MLD had vcrominent minima associated with both the Cyevems 
Nt. ¢ northeast dnd the anucycione in tie coupes. 
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IV. DEVELOPMENT OF AN AIRBORNE OCEAN SAMPLING 
STRATEGY 


There are two distinct considerations that are related to the successful areal 
mapping of oceanic variables: 1) The sample density in the survey region; and 2) The 


application of the objective analysis technique to map these samples. 


A. SAMPLING CONSIDERATIONS 

Given an ocean region to be synoptically mapped in some variabie, and the 
aircraft asset(s) and air-depioyabie sensors to measure this variable, what airborne 
sampling strategy should be empioved to acquire a representative sample of the field 
that will optimize the information returned with <he minimum amount of time and 
etfort expended to accurately map this variable? 

The ocean is a vast. temporally and spaually varying inhomogeneous domain. 
Following a simuiar development by Keilev (1976), this domain can be described by 


cn 


basic probabiity theory as a sample space, 5457.7 ee 


S,), of n observations where 


s = ¢ m = : 5 : bess * : , P 5 : » = ps ’ 2 * , * x ;| 2 x 9 
ig es Sy eres cess 1, 5 an ciement Of ene Sample space Ceimea OV TAUimGiee 


- 


longitude, and time. Each element can be further described by an area Vij 
(Vi (2), V5 (Zr Vi (2), j= 1,...,.m, of specific oceanic variables (biological, chemical, 
dynamic, or physical) varying with depth (z). Each V; j(2) can be defined as a 
continuous random variable which can have more than one potential value as the 
chance result of the sampling and that it is assumed to be representative of the oceanic 
variable being sampled. Due to the inherent capabilities of sampling instruments and 
systems such as the ADDAS, these oceanic variables are sampled at a specific depth 
increment and can be represented as discrete random variables. This provides a 
framework for the development of an “optimum” sampling strategy that can be applied 
in survey operations. 

Deciding which variable to sample is a rather simple task due to the availability 
of air-deployable sensors. The measurement of ocean temperature is most practical 
through the use of the SSQ-36 AXBT as it is a proven sensor in the Navy inventory. 
Synoptic mappings of the temperature field at selected depths can be converted to a 
sound speed field, using a climatological or representative salinity profile, for use in 
acoustic forecasting for ASW operations. Soon-to-be-operational AXCTD’s and 


currently available AXSV’s could also be used to map the sound speed field . 
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Having chosen the variable to be sampled, the “external” considerations of the 
spatial and temporal scales of the variable, and the “internal” considerations relating to 
aircraft, aircrew, and sensor limitations will be addressed. 

1. External Considerations 

a. Spatial Scales 

The ocean contains a multitude of spatial scales ranging from centimeters 
to thousands of kilometers. In sampling any oceanic variable, prior knowledge of the 
scales that are “significant” to a specific mission are required. Previous operations in 
the survey region, satellite imagerv, and/or climatology may be available. Additionally, 
a knowledge of the dynamical properties of the fieid must also be considered. 

The Rossby radius of deformation, A, is a length scale of fundamental 
importance in both atmospheric and oceanic dynamics for the behavior of rotating 
density-stratified fluids subject to gravitational restoring forces. For a homogeneous, 
shallow laver of fluid, it is the distance over which the gravitational tendency to render 
the free surtace flat is baianced by the tendency of the Coriolis acceieration to deform 
weemeutace (rediosky, 1979; p. #8). {t is the pertinent length scale for barotropic 
Phenomena {Emerv, er a/.. 1984). Since the ocean is a continuously stratified fluid, 
mermalsmode theory can te applied to a Stratified laver and each mode can pe 
represented by a Rossby radius. The barotropic or external Rossby radius (R) was 
defined above and is a function of only the water depth and latitude. The internal or 
nth mode baroclinic Rossby radius (Ri) depends on the vertical stratification and is a 
natural scale in the ocean often associated with fronts and eddies (Gill, 1982; p. 207). 
For mid-latitude oceans it is the order of 10 to 30 km. The horizontal dimensions of 
such quasi-geostrophic features should scale as the internal Rossby radius, where the 
scale 1s defined as the wavelength over 2% (Emery, et al., 1984). 

b. Temporal Scales 

Based on the spatial scales of the quasi-geostrophic features, the temporal 
scales of variability are not a factor in “synoptic” mappings. With phase speeds the 
order of 2 to 5 km/day, the distance that the features will have propagated over the 
duration of the survey flight is small compared to the sampling grid cell length. 
However, temporal scales may have to be considered if the survey domain is too large 


to be covered in a single flight. 
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c. Aliasing 
The well Known problem of aliasing occurring in time series analysis also 
occurs in spatial sampling. For a one-dimensional’ sampling, the Nyquist theorem 
States that some function can be reconstructed from samples taken not more than 
one-half wavelength apart. This can be applied to a two-dimensional field as well. 
Identify the highest frequency (shortest wavelength) of significant variability and 
sample at twice that frequency (one-half of that wavelength). If this is not feasible and 
a larger sampling distance is chosen, aliasing may be a probiem from scales that are 
smaller than the chosen sampiing wavelength and may account for a percentage of the 
noise in a sample value. 
2. Internal Considerations 
With the spatial scale of the features determined. and thus the maximum 
acceptable distance between observations to adequatelv represent the significant scales, 
the capabilities of the sampling platform and its sensors must be taken into account. 
These were presented as extrinsic limitations to the ADDAS in Chapter 2, and are 
applicable here. They are summarized: 
a. clircraft Limitations 
* Weight_and balance. Aircraft are constrained bv a maximum gross weight for 
take otf as prescribed bv the respective NATOPS Manual. There may have to be 
gee eso in the number of sensors to be deployed to remain within NATOPS 


e Airspeed and altitude. These ee affect fuel management. Airspeed and 
altitude will be governed by safety-of-flight and ADDAS/sénsor limitations. 


Additionally, physiological and safety factors must be considered. All Navy flight 
operations are governed by the OPNAVINST 3710.7 series (NATOPS General Flight 
and Operating Instructions). This specifies limits on mission endurance and crew rest 
policies. It may be further modified by individual squadron standard operating 
procedure (SOP). These will be the limiting factors in mission duration. 
b. Sensor Limitations 
¢ AXBT deployment time from the aircraft, a function of altitude. 


e AXBT fall duration through the water, 200 and 500 seconds for deep and 
shallow AXBI"s, respectively. 


e Three AXBT channels available. Two AXBT’s must not be deployed within 
line-of-sight transmission range of each other, a function of altitude. 
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B. GENERAL APPLICATION 

The thought process followed in developing an “optimum” sampling strategy 1s 
presented in flow diagram form (Figure 4.1). Each step is contingent on the previous 
step. Maximum observation spacing is generally given as a requirement to gain an 
adequate representation of the oceanic process. Aircraft maximum gross weight Is also 
a limiting factor. Trade-offs can easily be made between the size of the survey region 
and the maximum number of AXBT’s that can be carried. If the size of the survey 
region is then too large, two or more consecutive survey flights may be completed as 
Pemieecs tie tempor! scaies of the region are considered. One option. is to sampie the 
larger do.aain twice. This can be accomplished oy sampling in a “cneckerboard” 
pattern, with sampie spacing no more than twice the maximum observation spacing, 
with the first flight deploving AXBT’s in the “red” squares, and the second flight 


4 


depioying in the “dlack” squares. Another option is to compiete two or more surveys 
at a smailer-than-maximum spacing in adjacent subregions ot the sample domain. One 
drawback to sampiing with two or more ilignts is that the real-time aspect of che 
synoptic mapping will oe lost. This inay oe acceptable. however. depending on the 


fequlrements directing the survev. 


fe ah OBJECTIVE ANALYSIS TECHNIQUE 

The second area of consideration in developing an optimum sampling strategy is 
the proper application of the objective analysis technique. This will result in fields that 
are representative of the data acquired, and available in near-real-time with systems 
such as the ADDAS. 

Based on the Gauss-Markov theorem from statistical estimation and applied to 
atmospheric data by Gandin (1965), objective analysis was developed by Bretherton, ef 
al., (1976), with application to oceanic data. A more generalized implementation, as 
adduced by Carter and Robinson, (1981), is applied in the OPTOMA Program. 

1. Theory 

The basic assumption in objective analysis is that the field from which the 
observations were acquired has zero mean and a known covariance function. 
Following Bretherton, et al, (1976), and incorporating both a space and time 


variability, 0, is the value of a scalar variable @( x, y,1), at a general position p = 


p 
(x,y,t), from measurements @~, at N observations P; Ct— lee iN). Each 
measurement, @;, is presumed to consist of some true value, Q., and some random 


noise, E:: 
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Figure 4.1 
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9; = O(p; ) ot e3, (1 = ee, N ), 


where the errors, €, are uncorrelated with each other and with the true value, a but 


have a known variance, E: 


£0, = 0, §€ = £6;, for(i,j = 1,...N), 


and the overbar indicates an expected value (or average over a large number of 


realizations). The errors may be sampling error or instrumental error. Systematic or 


calibration errors are not permitted. 
The least squares optimum linear estimator of 0: 1S 
\ N N 
A= Yc.f V avg..) 
% = Cai! = aa sr a 


L= i j=l 


where 


+ E6.: 


— 


Ay = 0:0; = FCB - Bj 


i 


is the covariance matrix between all pairs of the sampled data, and A! is the inverse of 
A, and 


is the covariance between the estimated value 0, and the ith observation. The variance 


of the error in 0, is 


6 2 _ 1 
( p- 0, y= Cop : Pn vei Chi Aj, 


which requires only the location of the data points, a noise level E, and some 
correlation function F(R). 

Carter and Robinson, (1981), present the idea of applying the objective 
analysis equations at each interpolation point. Four specific parameters are identified 
that determine the selection of observations that will affect an interpolation. These are: 

e Maximum time increment between observations, 
Radius of influential points, R = [( Ax-c,Az > + Ay? ] 22 
e A dominant phase speed, c,, 


The maximum number of influential points. 


D1 


With application to synoptic mapping, only the radius of influence and number of 
influential points are used. It is best to use observation points near the interpolation 
point. By specifying a maximum range at which observations can effect an 
interpolation, this will eliminate the influence from distant points. Additionally, the 
specified maximum number of points limits only those observations within the radius 
that are the most highly correlated (closest to) an the interpolation point. The time 
increment and phase speed are disregarded because the propagation distance of the 
quasi-geostrophic features (phase speeds of the order 2 to 5 km/day), is negligibie over 
the duration of a survey flight (~8 hrs). The radius of influence is now onlyetne 
distance between the interpolation point and the observation R = (Ax7 + Ay" 2, 
2. ADDAS Software Application 
a. input Parameters 
Having acquired the observations. the following parameters must ove 

specilied before performing the objective analysis: 

¢ Some form of the correlation function. F(R ), 

o™ lhe percentage of the observed variance altiibuled todos Case 

¢ The maximum radius of influence. R/oy, 
* The umuting number of intluential observations, Vy, 
e The objective analysis interpolation points, Y, = {*p¥p)},(n =1.5/), 
e A maximum autocorrelation value, Anan 
These parameters may be incorporated into the program as constants based on 
previous experience or determined for each observation field. A flow diagram (Figure 
4.2) summarizes the decision process. These parameters are presented below. 

(1) Correlation Function and Noise Variance Level. 

There are two distinct uses for the correlation function in the objective 
analysis technique: 1) to determine the correlation between an interpolation point and 
an observation, and 2) to compute the elements of the autocorrelation matrix between 
all pairs of observations selected to influence an interpolated value. 

The correlation function can be defined in any number of ways 
depending on the user’s requirements. For a space-time correlation of asynoptic data, 


as is used in the statistical model, the form of the correlation function 1s 


F(R2) = ¥(1-@R2)exp( -BR2), (4.1) 


where 
2 ea )y? + (Ay- cyAt )*, the influential radius, squared, 
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Figure 4.2. Flow Diagram of Objective Analysis Inputs. 
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Cyr Cy! X and y phase velocity components, respectively, 


Ax, Ay: spatial difference between interpolation and observation points, 
At: time difference between observation and analysis times, 


a, B, ¥: parameters computed by an iterative nonlinear least squares fit to 
the function. 
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A simple, yet effective form is emploved in the ADDAS objective 


analysis soitware tor synoptic mappings. 
F( RO) = (1 =R2 lexo( = (4.2) 
where 


RS =0 Sone Roe PRT ~, the influential radius. squared. 


maxX? 


An advantage in using Equation 4.1 in synoptic mapping (let c, and C,, 


xX 
both equal zero) is that an estimate of che noise variance in the oGbservanons came 


made if not already known. First. a bilinear trend 1s removed (if desired) frommeaam 


observation. it is of the torm @:; = Ax, + By. +C, where A. 3B, and C are .emes 


i 4 
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specified or computed bv a ieast squares method. The mean, 4, and the variancejess 
of the observations are computed next, and the mean is removed from each 
observation to satisfy the basic assumptions. 


The distances between each observation and every other observation, 
I * 
Di = [04 - Meat vie, P17, i= 1a 


are computed and binned as if isotropic. Bin size should be greater than the minimum 
inter-observation distance ( Di;, 1*j) in order to accumulate a statistically significant 
number of values in each bin. The first bin should be from 0 to 1 km, to contain those 
Di when 1 = j, and other closely-spaced observations ( Dj; = 1km). The numbenen 
bins, mbin, should cover all the inter-observation distance scales. As each De is located 


in its respective bin, b, the total number of elements in that bin, NB}, 1s incremeniall 
by 1. The autocovariance is computed and summed for that bin ( CO, = » 59; ). 
Once all the pairs of observations have been binned, the final autocovariance for each 
bin, C, is computed by normalizing each C0, by the quantity NByo? 


(Cy, = COp/(NB,97)). 
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The midpoints of the bins, BM), for b = 2,...,nbin, are computed. 
These values, along with the total number of observations ( NB, ), the normalized 
correlations { C, ), the variance (o7), and the limiting number of observations ( N/2), 
are applied in an iterative nonlinear least squares fit to the correlation function 
(Equation 4.1). The values of a, B, y, and the Zero-crossing distance are computed. 
Some quality checks are made on the data. If the total number of values in either bin 
2, 3, or 4 is less than N/2, there are not a statistically significant number of 
observations used in computing the correlation, and the least squares fit !s not 
attempted. Additionally, the application of the least squares routine for this form of 
the correlation function requires at least three points to fit a curve, and these points 
must also occur before the second zero-crossing of the correlation values. The 95% 
confidence limits. based on a test for white noise and assuming a normai distribution, 
can be computed for each bin using the formula Cy = 1.96/( VB, 72) (Jenkins and 
Watts, 1968: pp. 187-188). 

If a successful fit is achieved. the signal and noise parts of the variance 
can be estimated by extrapolating the fitted correlation function to zero iag. Assuming 
perfect correlation at zero lag (C = 1.0) and that the variance contains some 


Meneentave Of 2O0ISE, if the extrapolated value, C..., is less than 1.0, the signal variance 


4 
can be computed, 


and the noise variance 1s 
jig ; 2 
a eee le On 


The noise variance must be specified as a fraction of the total variance for input into 
the objective analysis routine. 
(2) Maximum Radius of Influence. 

Ri yay 18 selected to be less than or equal to the zero-crossing distance 
of the correlation function. This ensures that only observations that are positively 
correlated with the interpolation point are used. 

(3) Maximum Number of Influential Observations. 

In surveys with a large number of observations, Nj, limits the 
number of observations within the radius that will be used in an interpolation. These 
observations will be those that are the most highly correlated with (i.e. closest to) the 


interpolation point. 
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(4) Grid Selection. 

For synoptic mappings in the OPTOMA Program, the survey domains 
are rectangular in shape, with observations distributed regularly over the domains. The 
objective analysis interpolation points, Y¥, = {x,,yp,}, (n =1,...,/), are organized in a 
rectangular pattern encompassing the domain, with equal spacing in both the local x 
and y directions. The grid spacing should be less than one-half of the zero-crossing 
distance of the correlation function to adequately resolve the scales of the sampled 
field. The interpolation points should lie within the dimensions of the domain tn order 
to provide enough observations for interpolation at the boundaries. 

(5S) Maximum Autocorrelation Vaiue. 

Amax 18 used in computing the autocorrelation matrix, A, at each 

interpolation point. Once the observations are selected that will influence an 


interpolation, the autocorrelation of each point with every other point, A::, is 


ij’ 
computed using the autocorreiation function F{ R) previously determined. If forga 


given j, the absolute value of A:, is greater than doy 


observation will be discarded. This removes observations that are too closely spaced to 


for any i, then them 


one another and heips maintain a diagonallv-deminant matrix. 
0. ine Objective Anatysis 
Once the input parameters are specified, the following steps are performed 
for each interpolation point, 
e Selection of the influential observations, 
e Generation of the autocorrelation matrix, 
¢ Computation of the interpolated value and error limit. 
These are discussed below. 
(1) Selection of Influential Observations. 
The influential radius is computed between the interpolation point and 


each observation. If the radius is less than or equal to RI,,,, then the correlation 
3 


x 
between the observation and the interpolation point C,; = F(R), is computed. Once 
the observations have been selected, they are sorted in order of decreasing correlation 
(increasing distance from an interpolation point). Only the first Nj;,, of the sorted 
observations and their respective correlations will be used. Next, the autocorrelations, 
Aji, 
and distributed observations with the interpolation point remain. 


between each of the selected observations are computed. Only the best correlated 
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(2) The Autocorrelation Matrix. 

The autocorrelations, AG, between the remaining selected observations 
are calculated using F( R), and become the elements of the autocorrelation matrix. 
The percentage of the variance due to noise is then added to each element of the main 
diagonal of A. This biases the the matrix towards diagonal dominance. The matrix is 
then inverted. This provides the weights to compute the interpolated value. The 
observation correlation vector C,,, between each selected observation and the 
interpolation point 1s computed using F7{ R ). 

(3) The {nterpolated Vaiue and Interpolation Error. 
The vaiues of the inverted autocorrelation matrix, A™'. and the 


correlation vector, C,,; are used to compute the interpolated value, B.,, 


NH. 
NV ay 
De wer A 


a i al ~ | a 
mo aa a 
and the RMS interpolation error 1s 
NV VY 
ee eee eC , AG; 


i= i 


. 


The interpolated value must now be be brought back to its original 
scale. If a bi-linear trend was removed, it must be added back, as well as the mean. 


The final result, R, 1S 
ae Oat ey tC; 


which is stored, along with En. for later contour plotting. The process is again 
executed for the remaining interpolation points. 
3. Possible Problem Areas 

From the numerous runs of the objective analysis that have been made in 
gaining an understanding of the technique and processing data, two interrelated factors 
have been found that determine the caliber of the objectively analyzed fields: data 
density and quality, and the noise variance estimate. 

Generally, the more robust the sample field, the better the quality of the 
resulting objective analysis. However, this may not always be the case. Too many 
selected observations too close to each other (high autocorrelation values) may result 


in an autocorrelation matrix that is not diagonally dominant, and an ill-conditioned 


i 


(singular or near-singular) matrix mav result. If this happens to be the case, the 
interpolated value may be grossly non-representative of the data it was computed from, 


and the RMS interpolation error will be high. By specifying 4 to be less than a 


max 
certain value (limiting the minimum distance between observations), this will limit the 
magnitude of the off-diagonal elements and result in a properly conditioned matrix, 
which should produce acceptable results. 

The estimated noise variance is another factor that can affect an interpolation, 
although the effect will not be as great as improperly specifying the maximum 
autocorreiation value. [fF the noise variance is specified to oe too small, especialigaam 
sparse data fields with a nigh noise variance, this mav also cause the autocorrelation 
matrix to oecome til-conditioned. As the noise is added to the main diagonal of the 
matrix. 1t also biases it toward diagonal dominance so a successiul inversion can 5é 
accomplisned. With a higher density of observations. smaller bin sizes can be used to 
compute the autoccrrelation vaiues and more points wiil be available for the least 
squares fit. This coutd result in a oetter estimate of that part of the varia 
attributabie to noise. 

Another name for the objective anaivsis technique ts optimal interpolation. 
With a xnowledge of the statistics ct the datavand a correlation function, the tracaem 
of the variance due to noise can be caiculated. These factors, along with a specified 
maximum autocorrelation value and the observations will result in an optimal 
interpolation of the observations to each interpolation point, and an acceptable 
representation of the sampled field. Failure to properly specifv these parameters may 
result in erroneous mappings due to the problems of non-singular matrices and 
observations that are too closely spaced. The effectiveness of a system such as the 
ADDAS for synoptic mappings is directly related to quality of the acquired data and 


the proper processing of this data. 


D. THE OPTOMA PROGRAM SAMPLING STRATEGY 
Applying the considerations discussed in the previous sections, and with the 
knowledge gained from earlier OPTOMA Program AXBT survey flights, the following 
key factors were considered in the development of the airborne sampling strategy in the 
NOCAL and CENCAL domains of the OPTOMA Program: 
e Size of region: The 250-km square covers the _maximum_area that can, be 


surveyed by a ey P-3 aircraft based at Moffett Field Naval Air Station, 
California, on a single mission. 
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e Number of AXBT’s: 96 AXBT’s is the maximum number that can be carried by 
the P-3 to remain within specified gross weight limitations and still have adequate 
fuel, including reserve, to complete a mission in this size region and at this 
distance from base, and to provide an adequate resolution of the mesoscale 
processes. 


e AXBT spacing: This is a function of many variables. AXBT’s are deployed at 
five minute. time intervals i 28 km spacing). . This increment has been 
determined from previous flights to Poe an Optimal coverage of the region, 
and it takes into account the spatial scale of oceanic features to be observed, 
ADDAS response, AXBT fall time from the aircraft, AXBT transmission time, 
and AXBT-aircraft line-of-sight transmission ae Or aircraft on-station 
groundspeed of 180 to 200 kts and altitude of 500 to [000 ft AGL. 
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V. SIMULATED APPLICATION OF THE SAMPLING STRATEGY 


The sampling strategy was applied to “simulated” AXBT survey flights in 
“synthetic” oceanic domains. The oceanic variability in these domains was synthesized 
from the superposition of linear planetary Rossby waves (LPRW’s). The simulated 
survey domains were 256-km by 256-km square, representing a nominal areal survey 
domain as in the OPTOMA Program. Synoptic mappings were generated from the 
simulated AXBT observations using the objective analysis technique. Some 
assumptions were made in applying the sampling strategy; these were: 

* No limit to aircraft endurance, airspeed, or payload capability, 
wi niimmited number of Ao s: 
2 No AXBT fatiures. 
* No limitations on the ADDAS, 
P iemporal scales or the elds were not cons:derew: 
Although not fully realistic, these assumptions provide a basis to exercise the oojective 


inalysis tecnnique on controlled data. 


A. DEVELOPING THE SYNTHETIC FIELDS 

A field of eddy-like features can be generated from the superposition of LPRW’s. 
Starting with the non-linear, hydrostatic, Boussinesq equations, and using a 
non-dimensionalization and perturbation expansion in terms of the Rossby number, a 


quasi-geostrophic potential vorticity (QGPV) equation can be derived: 


where 


sd = 
Dt dat Ox oy 


(The reader is directed to Gill, 1982, Ch.l!, for the complete Gerivation: mee 
governing equation (GE) for LPRW’s can be obtained by linearizing the QGPV 
equation, 

Op 


Civ Peace “Poy + B = 0 
at a Nee a 
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where, 
V7 () = GC /Ox? + 8° dy”. 
The assumed plane-wave soiution for the GE, 
Py = P(z) exp[i(ot-K(xcos® + ys) ) |, 

mere the real part, 

Dy = P(z) cos ( ot - kx - ly }, 
with 

A = « cos 9, and/ = x sin 6, 


is substituted into the GE and simplified, resuiting in 


ae — XpereeeS ae o oP{Z) 
pilot eV eein 2) | (ko 1) Ki G)| == —( — =)! 


OZ N~ 0O7 








which is the Ge for the vertical structure. Letting 


gives the resulting form of the equation, 


Oo 1 oP 


La 
dz Naz? BP? 


which is a Sturm-Liouville problem. Applying the rigid (kinematic) boundary 
condition, 0p/dz = 0, at z = 0, and at z = -D (where D is the water depth), there will 
be an infinite number of eigenvalues and eigenfunctions, where MW, are the eigenvalues. 
Solving for 6, the dispersion relation (DR) for the nth vertical mode is 

c.= aay 


—5——_2——., (n= 0,1, 2,...). 
ema) 0/229) 


The nth mode Rossby radius of deformation can be defined as A, = wo. In 
. e e e e 1 
dimensional form, the barotropic Rossby radius is hy = ¢y/ fp, where cy = (gD ) 44 : 
(wave speed) and fj = 2&2 sin(@) (coriolis parameter, and @ is the latitude). 
Assuming a constant Brunt-Vaisala (N*) profile, the baroclinic Rossby radii can be 


defined as ha = ND/ntjf,. By convention, 6, must be greater than or equal to zero; 


61 


therefore, cos ® must be less than or equal to zero, (n/2 5 9 S 3n/2), giving a 


westward, with a northward or a southward component, direction of propagation. 


B. THE SYNTHETIC Yizees 
The synthetic fields were calculated by the following general equation for WV 
different waves, each with N normal modes. Each wave is a suite of progressive modes, 


all having the same (x,/), and instantaneously in phase at t=0. 


iM V 
F(x. y,z.t)= >) » AQ, , @B(z) cos( G_, nt Rt x COsG em asm 
m= 1 ne 
where 
F: ts the value of the syntheric eid amseimes 7 2 
ee iS AN ardilrary amplitude function, 
bisa is a depth function, cos(%z/D), 
Jp: 1s the frequency lor the nrk normal mode, 
K_.: is the total wave number. -7. L. and 
Jo: is the propagation direction. 
The vnysical constants used in che caiculations were g = 9.8ms*, D = 2km. 
(eco and N= (0-< 10:85" giving fy = 9.07 = 1070 
B = 1.79 x 107!! (ms)7!, which are representative of a mid-latitude oceanic region as 


in the OPTOMA domains. The limits of the synthetic fields are -32 S x, y [288 km, 
to provide flexibility in sampling at the boundaries of the subdomain. The simulated 


sampling subdomain 1s within the region bounded by 0 S x,y S 256. 


€.. | PROCEDURE 
1. Synthetic Field Generation 

Two synthetic fields were generated through use of the general equation 
presented above. Each field was computed at t = 0 (Sn nt = 0) and z=0 
(ADE Ze 1s wach Amn = | for all waves and all modes. A summary of the specific 
parameters for each field is presented in Table 3. The origin of each wave contributing 
to the field is referenced to the origin of the subdomain, (0,0). “SYNFLD1” (Figure 
5.1) produced a standing pattern of eddy-like features with a wavelength of ~216 km 
(eddy diameter ~108 km). This field was scaled from 10 to 18 °C, representing a 
nominal SST range as observed in the OPTOMA domains. “SYNFLD2”’ (Figure 5.2) 
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Figure 5.1 Synthetic field, SYNFLDI and sample subdomain (dashed line), 
fepiesemune oof im ©. Contour interval 0.5°C. 
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Figure 5.2 Synthetic field 
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SYNFLD2 and sample subdomain (dashed line), 
representing MLD in meters. Contour interval 2.5 m. 
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was scaled to represent a mapping of mixed layer depth, with representative values 
from 20 to 70 meters. Again this was based on a similar range of MLD’s observed in 
the OPTOMA domains. 





TADPES 
See Corte rsPECIRI@*PARAMETERS 


Wield: SYNELD I Number of Waves: 2 Number of Yiodes: 2 

| Wave Origin Wavelength (km) eee) 

| mM = 7 m 
] (-42, 3) 150.0 i353 | 
. ter. 106) io0.0 26 | 

| fied? SYNPLD2 Number of “Yaves: 35 Number of Mlodes: 3 

| Wave.., Origin Waveiength (km) Cae! 

| | 2) 200.0 090 

| 2 432.-70 } iS.) L359 

| 3 iene! 150.0 180 
3 “99, £0}. Li 2.V ino 

| 5 (74, i110) *50.0 a0 


2. Simulated AXBT Surveys 

The simulated survey flights are based on the flight plan employed on 
OPTOMA20P, with flight tracks oriented parallel to the local x axis. 

a. SYNFLDI 
Progressions of six AXBT survey flights, each with a different AXBT spacing were 
conducted in SYNFLD1. Each different spacing was based on the reference spacing of 
Ax = Ay =28 km, as is used in OPTOMA Program survey flights. Four symmetric 
grids and two “checkerboard’ patterns were completed in each progression. 

(1) Progression 1. 

In this progression, four surveys simulated a symmetric x,y 
observation pattern with Ax = Ay = 14, 21, 28,and42km. These surveys were 
designated S1-14, S1-21, S1-28, and S1-42, respectively. The two remaining surveys 
were variations on the 21 and 28 km observations. Every other observation was used 
so that in each row, the samples were offset from each other, as the black squares on a 
checkerboard. These were designated $1-42X and S1-56X. 
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(2) Progression 2. 

In the second progression the same data were used, but a random 
position error, varying between = 4 km (~2 NM) was introduced to represent aircraft 
navigational error in AXBT locations. These represent one particular realization of the 
random position errors (used in the computations), and were designated 
S1-14N,...,S1-S56XN. 

bD.TSYNELDZ 

One survey was completed in this field wih Ax = Ay = 28 km. anda = 4 
xm position error. Again this represents a nominal OPTOMA survey flight. A 
progression of six Objective analysis computations were completed on the survey data. 
The two versions of the correlation function as presented in Chapter 1V¥ werouteces 
The input parameters. fraction noise variance and maximum autocorreiatiom values 
were varied tor €ach computation. 
3. Apptcation of Ghiective Anarysis 

rOr each simulated AXBT survey in a progression. the ODIlECtiVe anal 
technique was appiied as Chapter [V. itis briefly summarized: 

# Tne autocovariance function and the fitted correlation function are compugems 
and coniidence fimits at each point on the curve, pased on the 73% Coniid@eums 
interval {Or 4 nOnmal distribution are calewlarea: 

® From the fitted function, the fraction ot the variance due to noise is estimated 
and the zerc-crossing distance is determined. The noise variance, radius of 
influential points, interpolation point spacing, number of influential observations, 
and maximum autocorrelation value are then input into the objective analysis 
program. 

e The objective analysis technique is applied at each interpolation point, and 
quality checks on the value of the autocorrelation matrix determinant and 
maximum _ autocorrelation values are computed. The objectively analyzed fields 
and their RMS interpolation errors were computed and then mapped. 

4. Comparisons 

For each progression, the correlation (CORR), root mean square error 
(RMS), both the systematic (RMSS) and unsystematic (RMSU) root mean square 
error, and mean absolute error (MAE) were computed between interpolation points of 
each simulated survey and the true field (synthetic field observed at the interpolation 
points). Also each survey was compared with the _ reference survey 
(Ax = Ay = 28 km). The RMSS error is a measure of the linear bias in the objective 
analysis, while the RMSU error is a measure of the precision of the interpolated field. 
For a complete description of the formulas applied, the reader is directed to Willmott, 


et al., (1985). 
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D. SUMMARY OF RESULTS 
1. SYNFLD1 
The autocorrelation computations (Table 4) show an increasing estimated 
noise variance with an increase in sample spacing. The worst case was the 42 km 
S1-42 and S1-42N, for both progressions 1 and 2. The 


zero-crossing is also much shorter. This could be the result of the larger bin size used 


symmetric surveys, 


to compute the autocovariance values as the sample density was more sparse when 
compared to the other survevs. The relative magnitudes of the autocorrelation values 
in each progression are approximately the same. The observation positions, covariance 
42, S1-42X%, 
De 


i eo25): 


function, and fitted correlation function are presented for surveys Si-- 
$1-56X (Figures 5.5, 5.4, 5.5) and $1-42N, S1-42XN, S1-S6XN (Figures 5.6, 


TABLE 4 
AUTOCOVARIANCE AND CORRELATION FUNCTION 
| CALCULATIONS, SYNPFLDI 
| 
| Progression | 
| Survey Bin Sout ce Estimated Praction Feo 
SAE y signal Ses Noise SN: 
m aria e ariance istance 
S1-14 30 3.604 0.316 0.08] So9 
S1-21 30 3.455 0.483 0.123 57> 
$1-28 30 3.081 0.858 0.218 52.0 
S1-42 45 1.836 Ae 0.537 39.4 
$1-42X 30 B055 0.878 2723 So 
$1-56X 40 2.086 SoS 0.471 56.0 
Progression 2 
Survey Bin Estimated Estimated Fraction Zero 
Size igna nee Noise Crossing 
km) Variance ee Distance 
C (km) 
S1-14N 30 Be55 1 0.370 0.094 56.2 
S1-21N 30 3.584 0.343 0.087 553 
S1-28N 30 Se262 0.675 Oui 54.2 
S1-42N 45 DAs 1.852 0.467 43.6 
SI1-42XN 30 poo 0.343 0.087 56.4 
S1-56XN 40 Coe 1.011 0.256 60.5 
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Figure 5.5 SYNFLDI, Progression 1, Sample S1-56X. (a) Observation positions. 
(b) Autocorrelations (solid) and fitted function (dashed). 
95 % confidence limits (shaded). 
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Figure 5.7 SYNFLDI1, Progression 2, Sample_S1-42XN. (a) Observation positions. 
(b) Autocorrelations (solid) and fitted function (dashed). 
95 % confidence limits (shaded). 
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All objective analyses for SYNFLDI1 were interpolated to a regular 13 by 13 
grid with a spacing of 20 km. This spacing was chosen as it was less than half of the 
zero crossing distance of the reference survey, S1-28N. Input parameters were applied 
differently for each progression (Table 5). An estimated fraction noise variance of 9.05 
and a maximum autocorrelation value of 1.0 were used in Progression 1 to illustrate 
the results of improper specification of the objective analysis input parameters. The 
radius of influence was fixed at 50 km and the zero-crossing distance was 55 km. For 
Progression 2, the radius of tnfluence for each survey was the computed zero-crossing 
distance and the noise variance was as estimatec. The maximum autocorrelation value 
was selected to be 0.7, which iimits the distance between selected oObscrValiGmemes 
greater than 25 km. Equation 4.2 was the correlation function selectea@ =igmmae 
emploved, as in the ADDAS software. 

The objective analysis RMS interpolation errors are slightly larceraimm 
Progression 2 than in Progression !. but the range of vaiues for the determiniimee 
decreased bv three orders of magnitude (Taoie 35). Unis is a direct result Ciaiee 
specification of 4 smaiier maximum autocorreiation value ior Progression 2. ‘he 
appearance of the opjectively anaivzed fields remains consistent between (he wims 
orogressions :hrough surveys 31-28. Progression i. and 51-23N. Progression ime 
fields become more erratic with the larger observation spacing, $1-42, S$1-42X, and 
S$1-56X in Progression 1 (Figures 5.9, 5.10, 5.11), while maintaining a consistent 
representation of the field in Progression 2 (Corresponding surveys $1-42N, S1-42XN, 
and $1-56Xn, Figures 5.12, 5.13, and 5.14, respectively). Again, the cause of the erratic 
nature observed in Progression 1 appears to be the improper specification of the 
maximum autocorrelation value. The fraction noise variance may have a minor effect. 

The comparisons (Table 6) show that the correlation is very high for all 
samples in both progressions. However, samples $1-42N and S1-56XN in Progression 
2 show better correlation than S1-42 and S1-56X in Progression 1 when “proper” input 
parameters are specified. The magnitudes of the error parameters are generally 
consistent between the two progressions with the exception of the RMS error between 
$1-56X and S1-56XN. The systematic parts are nearly equal whereas the unsystematic 
part is almost five times greater for S1-56X. Since the RMSS is a measure of the 
precision of the interpolated field, and although Progression 2 (S1-56XN) contained 
errors in the positions of the observations, the proper specification, by the user, of the 


input parameters for the objective analysis appears to increase the precision. 
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HABLE 5 
OBJECTIVE ANALYSES SUMMARY, SYNFLD1 


Progression 1 





Input Parameters 


Sample Regs Vlaximum _ Maximum Fraction 
Number of Autocorrelation Noise 
tmaiwence Observations Value Variance | 
(Km) | 
ALL 50 7 0.7 0.05 
| 
Computations | 
Survey RMS Maximum Determinant 
Error Autocorreiation Maximum Viinimum 
{ fy } 
S1-14 0.01 0.918 O.185E-4 Delo 2-4 
Seal 0.01-9.03 Ona 0.590E-4 ). 2948-4 
| $1-28 U.U220815 0.344 ).987E-1 CS) 
51-42 O.01-0.29 Oey 0.850 0.9066 -i 
$1-42X 0.02-0.05 0.766 0.650E-1 O.874E-2 
51-50X% ).01-9,.22 ).344 rood O.514E-2 


Progression 2 


Input Parameters 


Sample Radius Maximum Maximum, Fraction 
of Number of Autocorrelation Noise 
Influence Observations Value Variance 
(km) 
S1-14N 56.15 ii 0.7 0.094 
S1-21N 55.09 ef 0.7 0.087 
$1-28N Ses 7 0.7 0.171 
$1-42N 43.56 Zi 0.7 0.467 
$1-42NX 56.39 a 0.7 0.087 
S1-S56NX 60.45 i 0.7 0.256 
Computations 
Surs ey RMS Maximum Determinant 
St Autocorrelation Maximum Minimum 
0 
14N ().04-0.07 0.700 0.107 0.109E-1 
2IN 0.04-0.04 0.696 0.418 0.115E-1 
28N 0.07-0.12 0.699 0.405 0.674E-] 
42N 0.31-0.49 0.446 ale Pale 6 
a0XN 0.04-0.12 0.697 0.174 0.130E-1] 
S6OXN 0.11-0.26 0.648 1.490 0.505 
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Figure 5.11. SYNFLD1, Progression,1, Sample 





1-56X. 


(a) Ovigeune analysis. Contour interval 0.5°C. 


) RMS interpolation error (%). 
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Figure 5.12 SYNFLD1, Progression 2, Sample $1-42N. 
(a) Objective analysis. Contour interval 0.5°C 
(b) RMS interpolation error (%). 


79 


200.0 


224.0 


tome) 


192.0 2240 ‘2560 


le SI-42XN, 
Oe OF 


aay 
%), 


Sam 


9 





160.0 


128.0 
Distance (km) 








9 
Nl 


—_ a 
— 
a 3, 
° ° eo 
9 
wy 
o ~ 
2 Be 


Sti a 2 

=A = i 

Sct ~ Sa 
Sete ~ _ - 


SEE) 
1490 
14.5 
I~ 
a 
/ 
32.0 64.0 96.9 











aa aie a = ; 

: ° ee — ZZN tN 

o9Se Oee OEl OOO! OB! 096 OF9 OFF 0 0952 OF? OZEI OO9t O821 096 O09 OZE OD 
(UIx) OUNISIY (uy) dour SIC 


160.0 


128.0 


Distance (km) 
ontour int 


96.0 
pieson x 
80 


32.0 64.0 
FDI Seco 

ective analysis. 

(b) RMS interpolation error ( 


0.0 
e 5.13 SYN 


(a) Ob] 


Figur 










WW 
e Sj) 
Ss 
14.5 


1920 2240 256.0 


w 
















160.0 





Ich 0 


Distance (km) 


96 0 






640 


“yy 
oe 


OU 


| Pgueece.. 6450 85960" [260 16010) 192.0 224:6 256.0 
Distance (km) 


Ng ey 
7 . ee =a a — 


206 U 


7 
V/ i a y | , : 


ww 
-— 


e2iQ 









192.0 






160.0 








Distance (km) 
1280 


96.0 








64.0 






13 7 
eaIN ee Wi 
0.0 32.0 64.0 96.0 128.0 160.0 192.0 224.0 256.0 
Distance (km) 


f2 











Figure 5.14 SYNFLDI, Progression 2, Sample S1-56XN. 
(a) Objective analysis. Contour interval 0.5°C. 
(6) RMS interpolation error (°%). 
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TABLE 6 


COMPARISON COMPUTATIONS, SYNFLDI 





Progression | 
Survey CORR MAE RMS RMSS RMSU 
| (°C) OC) CC) FC) | 
' 
S1-14 9,990 0.028 0,035 0,022 0.027 | 
Si-2! 0.998 0.104 Omen 9.040 0.124 | 
Sees 0.991 0.197 0,23} 0.090 267 | 
51-42 O65 0.360 U.518 0.038 Loy 
$1-=2X% G07 Los O.188 0.095 OnL02 | 
$1-36X% ioe Dees 1.019 On? 0.946 | 
Progression 2 
Survey ORR MAE RMS RMSS RMSU 
AS) ee ey oe 
S1-14N 9.996 O14} eo 0.060 ao? 
Si-2tN 0.998 OaelS al9 0.030 SO 
S1i-28N os ). 166 ea ).067 0.189 
$1-42N 0.994 ). 369 0.464 0.041 O.220 
SL-42XN 9.996 125 0.190 0.086 0.200 
Si-50XN 0.995 0.321 0.3383 0.530 ae 


2. SYNFLD2 

A single survey, representative of past OPTOMA survey flights, was 
completed. The autocovariance values and the fitted function were computed, resulting 
In a zero crossing distance of 74.9 km and a fraction noise variance of 0.162. The 
objective analysis interpolation points were selected to be at a 32 km spacing in a 
regular grid (Figure 5.15). 

Both versions of the correlation function, as presented in Chapter IV, were 
used with the objective analysis technique for comparison (Figure 5.16). The fitted 
correlation function (Equation 4.1) had a correlation of 0.838 at zero lag as computed 
by the least squares fit (corresponding fraction noise variance of 0.162). In computing 
the fitted function, it was assumed that there was perfect correlation (1.0) at zero lag in 
order to compute an estimate of the fractional noise variance. This value was then 
added to the main diagonal of the autocorrelation matrix. The ADDAS correlation 


function (Equation 4.2) had a correlation of 1.0 at zero lag. 
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Figure 5.16 Correlation functions as used with SYNFLD2. 
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In addition to different correlation functions, the specific input parameters 
that were varied in the objective analysis were the fraction noise variance and the 
maximum autocorrelation value (minimum distance between selected observations). 
These two parameters have been identified as problem areas in the objective analysis 
technique and required further investigation after the surveys on SYNFLDI. Two runs 
of three objective analysis computations were completed. The maximum 
autocorrelation value was changed for the two runs. It was selected from the graphs of 
the respective correlation functions (Figure 5.16), corresponding to distances of 25 km 
and 35 km. The two correiation functions were applied in the following manner. 

* The fitted function (Gauation 4.1) was used with the least squares caiculated 
values of the parameters @, §, 7. An estimated noise variance, as computed iv 
the the least squares fit, or 0.162 (Runs Al and Bl) was used. 

e The ADDAS function (Equation 4 .2) was used in two diiferent ways: |) with the 
fraction noise variance selected fo De 0. Q ot Ns ee and 82), and >) with the 
fractlon noise Variance as computed, 0.162 (Runs A3 and 235). | 

The zero crossing distance remained constant 1t 74.9 km and the maximum number of 
intluential points was chosen to be scven. Table * summarizes the inpul parameters 
mace results of the comparisons. 

All the objectiveiy analyzed fields show a smail range of RMS interpolation 
amors and <Spresent the true field well, with onty slight differences in the patterns and 
widely varying error reids. The smailest interpolation errors were observed on runs AZ 
and B2 (Figures 5.17 and 5.18, respectively), using the ADDAS function and zero noise 
variance for both the 25 km and 35 km minimum distances. The range of determinants 
for runs A2 and B2 was approximately one order of magnitude smaller than the other 
runs, however. 

In the comparisons, all the runs had a very high correlation with the true field, 
with runs A2 and B2 being the best correlated. The MAE, RMS and RMSU errors for 
runs A2 and B2 were one-half the value of the other runs, also indicating a better 
representation of the true field. Runs Al and B1 show the highest errors. 

3. Overall 

From the series of calculations presented, it appears that the higher the 
density of the observations, more observations will be available to influence an 
interpolation, and a better objectively analyzed field will result. However, the proper 
specification of the maximum autocorrelation value (minimum distance between 
selected observations) appears to be the key factor in the accuracy of the 
interpolations. Maximum autocorrelation values corresponding to one-third or 


one-half of the zero crossing distance appear sufficient for these data sets and should 
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CALCULATIONS, SYNPig@e2 





Objective Analysis Input Parameters 


Run Correlation Pracion Maximum. Corresponding 
Function Noise Autocorrelation Viinimum 
Variance Value Distance | 
! Al 4} 0.162 0.69 2) km 
| a oy 0.0 © 0.33 25 km 
| A3 a2 0.162 0.83 25 Km 
? Bl 4.1 9.162 0.50 3) Km 
B2 4.2 .Q 0. /0 2) kia | 
BS iz 0.162 0.70 So kaa 
Ohbiective Anaiysis Computations 
Run RMS Determinant. : 
EURO Maximum Minimum | 
(3 : 
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work well in regions of similar variability, using a similar correlation function. The 
noise variance does not appear to be a factor according to the SYNFLD?2 results, and 
could also be applied to the SYNFLD1 results as well. More work still needs to be 


done with the objective analysis technique and its application to synoptic mapping. 
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Figure 5.17 SYNFELD2, Run A2 (a) Objective analysis (m). 
(6) RMS interpolation error (%). 


87 





Ds 
O9S2 Ofce 


0 c6t 





0091 ORZI 096 
(UIy) Bue YsIq 





et ep ee 


26.0 


Distance (km) 


54.0 


S.0 


9.0 


O'O9T OBZ! 096 
(Wy) BOULISTG 





OF9 








: 


192.0 224.0 256.0 


64.0 96.0 1280 160.0 


32.0 


0.0 


Distance (km) 


Run B2 (a) Objective analysis (m). 


3 


D2 


MS interpolation error (%). 


NE 
R 


oa 
(b) 


Figure 5.18 


88 


VI. SUMMARY, RECOMMENDATIONS, CONCLUSIONS 


A. SUMMARY 

With the increasing array of aircraft deployable sensors, systems such as the 
ADDAS, ADAPS, and the Scripps System have proven invaluable in airborne ocean 
sampling and the real-time acquisition and analysis of oceanic data in the research 
community. In a matter of nours, aircrait surveys can cover an area that wouid 
require weeks with a research vessel. Synoptic mappings of the oceanic fields obtained 
from AXBT data can provide valuable information for the oceanographer and aiso 
have a direct application in ASW. 

Proper use Of such data acquisition systems requires a sampling strategy which 
acquires a representative sample ot the tield that will optimize the inrormation returned 
with the minimum amount or effort and time expended. A solid knowledge or the 
dynamics ot the oceanic variabies in the region must be considered, including both the 
f-@eret and Spatial scaies in order to sampie ihe field at a proper interval to 
adequately represent the significant features. Additionally, x<mowledge or the wrcratt 
Dlatiorm and its cupabilities, as weil as the limitations of the airborne sensors, must 
also be considered. With the data acquired, correct application of the objective 
analysis technique ts essential for the accurate representation of the synoptic mappings 


that will be generated. 


B. RECOMMENDATIONS 
1. A Proposed Ship-Aircraft Data Acquisition, Display, and Analysis System 
(SADADAS) 


The addition of real-time oceanic data acquisition and processing to a battle 
group can provide a stand-alone capability for the understanding of the oceanic 
environment and the tactical exploitation in ASW. Synoptic mappings in operation 
areas and at strategic transit points would give the warfare commander a “nowcast” of 
the oceanic temperature/sound-speed structure from which appropriate tactics can be 
employed. Furthermore, additional mappings at a later time in the same region can 
provide an assessment of the temporal variability in the region and also of the 


effectiveness of the tactics employed. 
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A system based on the ADDAS, ADAPS, and Scripps System is currently 
possible as these systems are research-proven, both in their hardware and software, and 
would not require a substantial lead time for development for fleet applications. The 
ADDAS 1s especially feasible as it uses totally off-the-shelf hardware. The ADDAS 
software is based on shipboard data acquisition software developed in the OPTOMA 
Program so a ship/aircraft software package is already available. The software suites 
of both the ADDAS and ADAPS are written in Hewlitt-Packard (HP) BASIC and 
could easily be transrerred to other HP-8ASIC systems. A software suite could be 
incorporated into the Geophysics Fieer Mission Program Library (GFMPL), and with 

the addition of a digitizing unit, could provide an interim real-time data acquisition and 
analysis capaoilitv until a comprehensive system 1s available. 

a. System Description 

The main hardware components of the SADADAS would be a2 
multi-channel digitizing unit and a multi-tasking computer. The digitizing unit snouid 
be easily selectable setween the array of air-deployable sensors (AXBTs. AXS ve 
AASVs, AXCTDs), as well as expendable sensors such as the XBT and <@C Sime 
a muitichannei system. the simultaneous processing of two sensors would require a 
computer “vith muiti-tasking capability. Both of these components are ‘veil within the 
available technology on the market today. Additionally, a printer and hard-copy 
graphics capability should be included. The HP-9020 and related hardware, currently 
being deployed with the Tactical Environmental Support System (TESS) in the Naval 
Oceanography Community, could provide the hardware capability with only the 
addition of a suitable digitizing unit. However, TESS uses FORTRAN-based software 
so the programs would have to be translated from HP-BASIC. 

The software suite would contain a menu-driven acquisition and storage 
program for each expendable sensor, both air and ship deployable. A program to 
generate the standard Navy and WMO AXBT- and XBTJJXX messages as the profiles 
are digitized should be included. This would result in less man-hours being used in 
manually selecting the inflection points from a hard-copy profile and provide more 
accurate profiles for archiving at Fleet Numerical Oceanography Center (FNOC). 
Another menu-driven software package would select which field is to be displayed. The 
data would be read from storage and objectively analyzed to produce the synoptic 
mappings. Hard-copy output would then be available to the warfare commander and 


oceanographer for tactical briefing and planning. 
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b. System Employment 
As the name SADADAS implies, the proposed system will be able to be 
deployed in land-based ASW aircraft as well as in ships with deployed air assets. These 
systems are presented below. 
(1) Land-Based Aircraft. 

The aircraft variant of the system would be employed in much the 
Same manner as in the OPTOMA Program aerial survevs as presented in Chapter ITI. 
The hardware and operator would be physically located onboard the aircraft so the 
system must be portable and compact. A system such as this could prove particularly 
attractive in regions of repeated coverage or continuous on-station operations. 
Although there will be some time lag between the initial survey flight and the 
availability of data to subsequent flight crews, synoptic mapvings of the oceanic 
Variability would be available for tactical planning and could be updated reguiariy cn 
following flights. A data base could be buiit up and archived at the local ASW 
Operations Center (ASWOC) for future missions. A system such as this would 
interface with the (TESS rather well. 

(2) Denioved Aircraft Assets. 

A shipboard variant of the ADDAS could oe developed with ttle 
effort and would provide another facet to real-time data acquisition. The hardware 
and operator would be stationed either in the carrier ASW module or in the Combat 
Information Center (CIC) in destroyers and battleships. AXBTs could be deployed by 
S-3 Viking carrier-based ASW aircraft and the SH-2 Sea Sprite, SH-3 Sea King, and 
SH-60 Seahawk ASW helicopters deployed in frigates and destroyers in the battle 
group. The AXBT audio signal would be transmitted to the ship by the aircraft-ship 
data link. AXBT locations would be provided by either the shipboard air controller or 
by the aircraft onboard navigation and transmitted to the ship. The AXBT audio 
signal would be received from the shipboard sonobuoy lofargram recorder and fed into 
the digitizing unit. Two or more aircraft could be utilized to perform the survey and 
provide larger areal coverage of the region. With the survey complete, near real-time 
synoptic mappings would be available for operational use and tactical decisions. 
Again, this capability could be accomplished with the TESS now being deployed. 

Factors such as the AXBT loadout, data link transmission range, 
navigational accuracy of AXBT positions, and speed and endurance of the helicopters 
need to be investigated more fully to determine the feasiblity of a shipboard data 


acquisition capability and its employment. 
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c. Operationally Applicable Synoptic Fields 
The following fields are presented as tactical decision aids for the warfare 
commander: 


¢ Sea Surface Temperature. This will indicate the surface manifestations of fronts 
and eddies that may have tactical significance in ASW. 


e Temperature/Sound speed at selected depths. This will determine the depth 
extent of the surface features and provide mappings of tactically significant 
oceanic structure, especially with respect to acoustic propagation. 

e Mixed layer depth. This will provide the warfare commander with the depth of 
the surface duct (if one Cu This has particular application in surface ship 
oow-mounied sonar, towed-array sonar, and sonobuoy employment. The 
acoustic properties of the mixed laver will ultimately determine the transmission 
ranges in the surface duct for bow sonars., A knowlédge of the mixed layer depth 
can provide the towed-array operator with an optimum depth to employ the 
array. Sonobuoys are limited by the depth of their hycrophones. The shallow 
sonoouovs have a hydrophone depth of 90 feet (~ 27 meters). With variability in 
the mixed layer, the hydrophone mav be in the layer at one location while below 
the layer at another. “This will affect the transmission and receiving capabilities 
of the sonobuoy. 

* Subsurface ducts and sound channels. This will have a similar application as in 
the depth of the surface duct (mixed layer}. Subsurface sound channeis may 
atfect the the acoustic propagation in the Survey region. 

* Vertical sections of temperature/sound speed. Vertical sections through a survey 
region can indicate temperature/sound speed gradients along a specific track from 
the surface to the depth or the survey (normally 300 m}. Sound speed pneu. 
could also be zenerated for input to acoustic prediction models to get the 
transmission loss and acoustic paths along the track. 


C. CONCLUSIONS 

With the capability of research-proven systems and the availablity of off-the-shelf 
hardware and expendable sensors, the development of a SADADAS for the operational 
Navy is currently feasible. The addition of a multichannel digitizing unit and data 
acquisition software to the TESS could provide the Naval Oceanographer and the 
warfare commander with a near-real-time, stand-alone ocean sampling and analysis 
system that can enhance tactical readiness. 

There is an emerging need for ocean sampling and synoptic mapping of the 
mesoscale variability in naval operations, especially in antisubmarine warfare. As 
submarine-radiated noise has been decreasing as advances in quieting technology are 
implemented, the warfare commander will have an ever increasingly difficult task of 
localizing and tracking submarines. Real-time synoptic mappings from AXBT surveys 
can provide the Naval Oceanographer and the warfare commander with a 
three-dimensional representation of the thermal (and hence acoustic) structure of the 
submarine’s environment. Standard ASW tactics can be modified and new ones 


developed as the fleet is educated in the variability of the oceans. 
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A most valuable, though unanticipated result of the ADDAS deployment in the 
OPTOMA Program, was that Navy aircrews became aware of the variability of the 
ocean. Usually only one or two AXBTs are deployed in an exercise or mission. These 
AXBTs may not be representative of the ocean structure in the region yet are used to 
compute acoustic parameters which are applied to the sensors. A system such as the 
SADADAS could provide aircrews and shipboard tacticians with a new perspective of 
the ocean, mesoscale variability and its effects on operations. It is time for such a 


system to ve developed further and evaluated in actual fleet operations. 
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APPENDIX A 
ADDAS HARDWARE COMPONENTS 


The following components are used in the current version of the ADDAS. Prices 
listed are the approximate retail cost as of July 1986. Svstem build-up requires 
approximately 40 hours at the GS-7 level. Software testing/maintenance requires 


about 10 hours at the GS-9 level. 


HP9816 MC68000 16-bit Microprocessor $5,000.00 
w/HP912! 3.5 inch microtloppy dual disk drive 
Sippican Vark 9 Front-end Digiuzing Unit 5,000.00 
Pes Wireuie. board “00.00 
HP7470A Two-oen Plotter 1,200.00 
Avionic Instruments Modei ZA125-iA Inverter 3,000.00 
Vil Stereo ii-7 1 yy ideocassetve seconde. 800.00 
TeK(rOnixes! b Oscilescome 750.00 
Portable Rack- Mount Storage/Transit Case 500.00 
ARNAY R40 LORAN-C Unit (optional) 400.00 
w/Antenna and cable 
Assorted electronic supplies 200.00 
Computer interface cables 
System wirin 
Assorted knobs and switches 
Total Cost $17,550.00 
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